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Abstract 



We study the variable bottom generalized Korteweg-de Vries (bKdV) equation dtu — —d x (d^u + f(u) — b(t, x)u) , 
where / is a nonlinearity and b is a small, bounded and slowly varying function related to the varying depth of a 
channel of water. Many variable coefficient KdV-type equations, including the variable coefficient, variable bottom 
KdV equation, can be rescaled into the bKdV. We study the long time behaviour of solutions with initial conditions 
close to a stable, 6 = solitary wave. We prove that for long time intervals, such solutions have the form of the 
solitary wave, whose centre and scale evolve according to a certain dynamical law involving the function b(t,x), plus 
an _H" 1 (R)-small fluctuation. 

1 Introduction 

We study the long time behaviour of solutions to a class of Korteweg-de Vries type equations, which we call 
the variable bottom generalized KdV equation (bKdV). These equations are of the form 



where / is a nonlinearity and b(t, x) is a real function. When f(x) = x 2 , the bKdV is related to an equation for 
the bottom of the channel appearing in the derivation of the KdV from shallow water wave theory. Examples 
of possible choices for the nonlinearity are f(u) = u 2 , the Kortweg-de Vries (KdV) from shallow water wave 
theory; f(u) — u 3 , the modified KdV (mKdV) from plasma physics; and f(u) = u p , the generalized power 
nonlinearity KdV (gKdV). When b = 0, JIJ reduces to the generalized Korteweg-de Vries equation (GKdV) 




(1) 



d t u = -d x (d 2 x u + f(u)). 



(2) 
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The KdV is obtain by unidirectionalizing the small amplitude, long wave/shallow water limit of the two 
dimensional water wave system with a constant bottom. The first such derivation was given by Korteweg 
and de Vries |32| over a century ago in an attempt to explain the existence of solitary waves of permanent 
form in a shallow channel. Numerous authors have improved the formal derivation using either asymptotic 
expansions or Hamiltonian methods ^|^]. Schneider and Wayne have given a rigorous proof of 
the validity of the KdV in approximating the water wave system in the KdV regime over time intervals of 
01. The KdV also appears in algebraic geometry. A nice survey of the KdV and it's relation to algebraic 
geometry is given by Arbarello P^. 

A remarkable property of the GKdV is the existence of spatially localized solitary (or travelling) waves, 
i.e. solutions of the form u = Q c (x — a — ct), where a£l and c in some interval /. When f(u) = u p and 
p > 2, solitary waves are explicitly computed to be 

Q c (x) = c~^Q{ci(x - ct)), 

where 

<3W = (£±i) A (co,h(^ 

It is generally believed that an arbitrary, say solution to equation © eventually breaks up into a 

collection of solitary waves and radiation. A discussion of this phenomenon for the generalized KdV appears 
in Bona [H]. For the general, but integrable, case see Deift and Zhou 19 . 

As channels with constant bottom do not exist, it is of interest to know how solutions initially close to 
a solitary wave behave as the wave propagates over channels with a variable bottom. Derivations of KdV- 
type equations when the bottom varies slowly have been presented by numerous authors. See, for example 
[lfil 1441 12~71 1331 14 7j . The resulting equations are the KdV with variable coefficients depending on the variable 
bottom. These derivations are non-rigorous and agree to the leading order in the bottom length, i.e. in 
order sup \d x h(x)\. We assume a depth h(x) of 01 with length scale lh- Consider solutions of the water wave 
system with wavelength scale l\, wave amplitude scale l n and fluid velocity scale l u . Then, if these scales are 
related as l\ = 0e^ , l v — 0e, l u — 0e and lh = o(e^), the leading equation for the wave amplitude, after an 
additional rescaling of the time variable, is (see equation (74) in |44p 



dtv = -r (2 + Jy + \d x 



Here T is the anti-symmetric operator r := | (c(x)d x +d x c(x)), with c(x) — y/ gh(x) (g is gravitational 
acceleration), and er](e~^x, e~it) is the surface elevation measured from the flat interface y — 0. We assume 
h = ho + hi with \hi\ << ho and ho a nonzero constant. Dropping terms of 0d x h in the above equation for 
r/, and changing variables as ^(x, i) — v{y, t), where y = x — and c$ — \fgh~o, leads to an equation for v: 



d t v = -dy -v + —cv 1 + -chfdiv 

\ e 4ft 6 y 



To order 0h±, solutions of this equation and solutions u{x,t) to the bKdV with nonlinearity f(u) = u 2 and 

b(x,t) = i (^=c§ho (x + - c ^j 
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are related by the transformation 

, 4h ( VE Ve \ 

v{y, t) = —u — — y, — — t . 
6c o \c>h 4 h J 

In a wider range of parameters one should add more complicated, in particular, nonlocal terms to ffl. We 
expect that the modified equation can still be treated by the methods developed in this paper. 

Similarly, in many other instances in mathematics and the sciences where the GKdV arises from an 
approximation of more complicated systems, the effects of higher order processes can often be collected into 
a term of the form b(t, x)u. Our main result stated at the end of the next section gives, for long time, an 
explicit, leading order description of a solution initially close to a solitary wave. 

We assume that the coefficient b and nonlinearity / are such that QJ has global solutions for H 1 (M) 
data and that with b — possesses solitary wave solutions. We discuss the latter assumption in Section 
13 Here we mention that the literature regarding well-posedness of the KdV (b — 0, f(u) — u 2 ) is extensive 
and well developed. Bona and Smith proved global wellposedness of the KdV in H 2 (R). See also [2"%] . 
More recently, Kenig, Ponce, and Vega [5U] have proved local wellposedness in if s (R) for s > — | and global 
wellposedness in iJ 1 (R) for s > 1. Similar results are available for the gKdV |29|. More recently, local 
wellposedness results in negative Sobolev spaces for the KdV have been extended to global wellposedness 
results. See ^)EJ- We are not aware of a wellposedness result for the bKdV in Hence, in the 

next section we give a global wellposedness result, whose proof (see Appendix EJ uses results of [211, an d 
perturbation and energy arguments. We conjecture that global wellposedness remains true for b £ C 1 
bounded and subcritical nonlinearities. 

Soliton solutions of the KdV equation are known to be orbitally stable. Although the linearized analysis 
of Jeffrey and Kakutani (20) suggested orbital stability, the first nonlinear stability result was given by 
Benjamin (2). He assumed smooth solutions and used Lyapunov stability and spectral theory to prove his 
results. Bona 0] later corrected and improved Benjamin's result to solutions in H 2 (R). Weinstein 05] used 
variational methods, avoiding the use of an explicit spectral representation, and extended the orbital stability 
result to the GKdV. More recently, Grillakis, Shatah, and Strauss [21] extended the Lyapunov method to 
abstract Hamiltonian systems with symmetry. Numerical simulations of the soliton dynamics for the KdV 
were performed Bona et al. See jUJ 111)1 |SJ . 

For nonlinear Schrodinger and Hartree equations, long-time dynamics of solitary waves were studied by 
Bronski and Gerrard Frohlich, Tsai and Yau [221, Keraani [3J, and Frohlich, Gustafson, Jonsson, and 
Sigal For related results and techniques for the NLS see also [T31 IT31 031 1571131)1 1331 WR ETH 155) . 

In our approach we use the fact that the bKdV is a (non-autonomous, if b depends on time) Hamiltonian 
system. As in the case of the nonlinear Schrodinger equation (see we construct a Hamiltonian reduction 
of this original, infinite dimensional dynamical system to a two dimensional dynamical system on a manifold 
of soliton configurations. The analysis of the general KdV immediately runs into the problem that the 
natural symplectic form oj is not defined on the tangent space of the soliton manifold. In the case of the 
mKdV (/(it) — u 3 ), the sympletic form is well defined on the tangent space because of the special structure 
of the solitary wave Q c , and hence Dejak and Jonsson |20] were able to prove long time dynamics of solitary 
waves in this special case. 

To address the problem regarding the symplectic form, we introduce a family of symplectic forms u) a 
parametrized by a small parameter a > 0, and approximating u>. We use the small parameter to control 
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the errors generated by this approximation. This approach works, except at one crucial step: the resulting 
lower (coercivity) bound on the Hessian of the energy functional is too weak to close our energy estimates. 
To remedy this we show that the weak bound comes from the directions in which we regularized lu; on the 
orthogonal complement the lower bound is sufficiently good. Hence, we decompose a general tangent vector 
into "bad" and "good" directions, and use precise information about the "bad" directions to considerably 
improve the upper bound (involving the nonlinearity) and close the energy estimates. 

In the next section we formulate our assumptions, state the main result and describe the organization 
of the paper. All L 2 (R) and Sobolev spaces used in this paper, except those in Section 0] are real. 

Acknowledgements 

We are grateful to J. Bona, R. Pego and Zhou Gang for useful discussions. 

2 Preliminaries, Assumptions, and the Main Result 

We begin with the following global wellposedness result proven in Appendix [S] See the appendix also for 
the definitions of the norms of b used in the following theorem. 

Theorem 1. Let uq 6 H 1 (R). For small enough 



there is a unique, global solution u G C(R, _ff 1 (M)) to Q with f = u 2 . With modification of the norm ||6||xt, 
the result continues to hold for f — u 3 and f — u 4 . 

The bKdV can be written in Hamiltonian form as 



where H' b is the L 2 (R) function corresponding to the Frechet derivative dH\, in the i 2 (R) pairing. Here the 
Hamiltonian iJf, is 



where the function F is the antiderivative of / with F (0) = 0. The operator d x is the anti-self-adjoint 
operator (symplectic operator) generating the Poisson bracket 



PWxt := \\b\ 



■4,oo 



+ INI^ i§? + ll&'IUi^. 



d t u = d x H' b (u), 



(3) 





defined for any F, G such that F', G' € (K). The corresponding symplectic form is 
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defined for any v±,V2 £ Here the operator d x 1 is defined as 

d~ x v{x) := / v(y)dy. 



Note that d x 1 • d x — I and, on the space {u £ L 2 (R) \ udx — 0}, d x 1 is formally anti-self-adjoint with 
inverse d x . Hence, if v\(x) dx = 0, then ui(vi, V2) — vi(x)d~ 1 V2(x) dx. 

Note that if 6 depends on time t, then equation @ is non-autonomous. It is, however, in the form of a 
conservation law, and hence the integral of the solution u is conserved provided u and its derivatives decay 
to zero at infinity: 

d f°° , 

u dx = 0. 



dt 

There are also conserved quantities associated to symmetries of with 6 = 0. The simplest such corresponds 
to time translation invariance and is the Hamiltonian itself. This is also true if b is non-zero but time 
independent. If the potential 6 = 0, then is also spatially translation invariant. Noether's theorem then 
implies that the flow preserves the momentum 

P{u) :=\\\ut L ,. 

In general, when 6 7^ the temporal and spatial translation symmetries are broken, and hence, the Hamil- 
tonian and momentum are no longer conserved. Instead, one has the relations 

(d t b)u 2 dx, (4) 



|p(») = ifW 



(5) 



where b'(t,x) := d x b(t,x). For later use, we also state the relation 



-- J bu 2 dx = J -u 2 d t b + b> (uf(u)--(d x u) 2 -F(u)) -b"ud x udx. (6) 



Assuming Q is well-posed in H 2 (M.), the above equalities are obtained after multiple integration by parts. 
Then, by density of H 2 (M.) in the equalities continue to hold for solutions in H 1 (M.). To avoid these 

technical details, we assume the Hamiltonian flow on i/ 1 (M) enjoys (JSJ and iJBJ. 

Consider the GKdV, i.e. equation (J2J. Under certain conditions on /, this equation has travelling wave 
solutions of the form Q c (x — ct), where Q c a positive H 2 (R) function. Substituting u = Q c (x — ct) into the 
GKdV gives the scalar field equation 

-dlQ c + cQ c - f(Q c ) = Q. (7) 

Existence of solutions to this equation has been studied by numerous authors. Sec 40, 3j. In particular, in 
[Sj, Berestyki and Lions give sufficient and necessary conditions for a positive and smooth solution Q c to 
exist. We assume g := — cu + f(u) satisfies the following conditions: 

1. g is locally Lipschitz and g(0) = 0, 
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2. x* := ini{x > | J* g(y) dy} exists with x* > and g(x*) > 0, and 

3. lim s ^ ^ < -m < 0. 

Then, as shown by Berestycki and Lions, (0) has a unique (modulo translations) solution Q c £ C 2 , which 
is positive, even (when centred at the origin), and with Q c , d x Q c , and <9 2 Q C exponentially decaying to zero 
at infinity (d x Q c < for x > 0). Furthermore, if / is C 2 , then the implicit function theorem implies that 
Q c is C 2 with respect to the parameter c on some interval Iq C M+. We assume that x m 9"Q c € i 1 (R) 
for n — 1,2,3, m = 0,1,2 and that J_ oo d c Q c dx ^ 0. The first assumptions are needed for continuity 
and differentiability with respect to c of integrals containing <9™<3 C , and the last assumption is made for 
convenience. When £™ dx = 0, unboundedness of d x does not present problems (see |2U|1. 

The solitary waves Q c are orbitally stable if S'(c) > 0, where 5(c) = P(Q C ). See Weinstein for 
historically the first proof for general nonlinearities. Moreover, in [21], Grillakis, Shatah and Strauss proved 
that 8'(c) > is a necessary and sufficient condition for Q c to be orbitally stable. In this paper, we assume 
that Q c is stable for all c in some compact interval I C Iq, or equivalently that S'(c) > on I. For f(u) = u p , 
we have S'(c) — 4 ^ p 1 ^ ||Q c= i|j L2 , which implies the well known stability criterion p < 5 corresponding to 
subcritical power nonlinearities. 

The scalar field equation for the solitary wave can be viewed as an Euler-Lagrange equation for the 
extremals of the Hamiltonian Hb = o subject to constant momentum P(u). Moreover, Q c is a stable solitary 
wave if and only if it is a minimizer of Hb=Q subject to constant momentum P. Thus, if c is the Lagrange 
multiplier associated to the momentum constraint, then Q c is an extremal of 

A ca (u) := H b=0 {u) + cP{u) (8) 

= J -{d x u) 2 + -cu 2 -F{u)dx, 

and hence A' ca (Q c ) = 0. 

The functional A CQ is translationally invariant. Therefore, Q ca {x) :— Q c (x — a) is also an extremal of A CQ , 
and Q c (x — ct — a) is a solitary wave solution of with 6 = 0. All such solutions form the two dimensional 
C°° manifold of solitary waves 

M s :={Q ca |ce/,a€M}, 
with tangent space TQ ca M s spanned by the vectors 

Cca : = daQca = ~-d x Q C a and Cca : = dcQca, (9) 

which we call the translation and normalization vectors. Notice that the two tangent vectors are orthogonal. 

In addition to the requirements on b that (JJJ be globally wellposed, we assume the potential b is bounded, 
twice differentiable, and small in the sense that 

\d?d™b\ < e a e?e™, (10) 

for n = 0, 1, m = 0, 1, 2, and n + m < 2. The positive constants e a , e x , and et are amplitude, length, and 
time scales of the function b. We assume all are less than or equal to one. 

Lastly, we make some explicit assumptions on the local nonlinearity /. We require the nonlinearity to 
be k times continuously differentiable, with /W bounded for some k > 3 and /(0) = /'(0) = 0. These 
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assumptions ensure the Hamiltonian is finite on the space 77 1 (]R) and, since Q c decays exponentially (see 
[3]), exponential decay of f(Q c ) and f'(Q c )- 

We are ready to state our main result. 

Theorem 2. Let the above assumptions hold and assume S'(c) > for all c in a compact set I C Iq. Let 
< s < |. Then, if e a < 1, e x < 1} £o << (ea^) 2,5 are smo/( enough, there is a positive constant C such 
that the solution to £JJ) with an initial condition uq satisfying wlQ ca eM s \[Uq — QcaWn 1 — e o can be written 
as 

u(x,t) = Q c (t){x - a{t)) + £(x,t), 

where \\£(t)\\ H i < &t s a e x for all times t < C(e% + e x + e^e^.) -1 . Moreover, during this time interval the 
parameters a(t) and c(t) satisfy the equations 

C-b(t,a) \ ,0 2s 2s, 2 

where c is assumed to lie in the compact set I. 

Sketch of Proof and Paper Organization. To realize the Hamiltonian reduction we decompose functions in a 
neighbourhood of the soliton manifold M s as 

U = Qca + £ 

with £ symplectically orthogonal to TQ ca M s , i.e. £-Ld~ 1 TQ ca M s . Unfortunately, since d~ 1 TQ i , a M s <2 L 2 (R), 
such a decomposition is ill-defined for £ 6 To overcome this difficulty we construct in Section 01 an 

approximate symplectic form 

1 f°° 

v{vi,v 2 ) = - / vi(x) (7C Qq - JC Qa *) v 2 (x)dx, 

^ J — oo 

where a > and K,Q a is a bounded operator regularizing the unbounded operator c?" 1 in certain directions. 
We show that there is an eo > such that if the solution u satisifes the estimate inf<g ca ||u — QcaWn 1 < e 0i 
then there are unique C 1 functions a(u) and c(u) such that u — Q c («)o(«) + £ with £ 6 (K.Q a TQ ca M s ) ± . 

With the knowledge that the symplectic decomposition exists, we substitute u = Q ca + £ into the bKdV 
(|T|l and split the resulting equation according to the decomposition 

L 2 {R) = JC Qa T Qca M s © (K.Q a TQ ca M s )' L 

to obtain equations for the parameters c and a, and an equation for the (infinite dimensional) fluctuation 
£. In Sections 0] and we establish spectral properties and an anisotropic lower bound of the Hessian A" a 
on the space (K.Q a TQ ca M s ) ± . Using these properties we orthogonally decompose £ into a "bad" part and 
a "good" part £ g , where is colinear with the minimizer -q of ||£||^ 2 (£,£q£). In Section|S]we isolate the 
leading order terms in the equations for a and c and estimate the remainder, including all terms containing 
£ = £f, + £ g . We use the special properties of the minimizer 77 to obtain a better estimate on the nonlinear 
terms containing £&. 
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The proof that H^Hjji is sufficiently small is the final ingredient in the proof of the main theorem. The 
remaining sections concentrate on proving this crucial result. We employ a Lyapunov method and in Section 
0we construct the Lyapunov function M c and prove an estimate on its time derivative. This estimate is later 
time maximized over an interval [0, T], and integrated to obtain an upper bound on M c involving the time 
T and the norms of £b and £ g . This anisotropic upper bound is considerably better than an isotropic bound. 
We combine this upper bound with the anisotropic lower bound M c > Ca W&Wh 1 + @ llCsll m > wrnc h follows 
from the results of Section[51 and obtain and inequality involving the norms of £b and £ ff . Here a is the small 
regularization parameter mentioned in the introduction and will be taken small, in fact a = (e a e x ) s . This 
inequality implies upper bounds on and provided 11^(0)11^1 is small enough, via the standard 

argument given Section[S] We substitute this bound into the bound appearing in the dynamical equation for 
a and c, and take e a e x and eo small enough so that all intermediate results hold to complete the proof. □ 



3 Modulation of Solutions 

As stated in the previous section, we begin the proof by decomposing the solution of Q into a modulated 
solitary wave and a fluctuation £ : 

u(x,t) = Q c (t) a {t){x) +£0M)> (11) 

with a, c, and £ fixed by an orthogonality condition, which we now describe. Ideally we would like to take £ 
orthogonal to JCTQ aa M Sl where JC is the symplectic operator defined on absolutely continuous functions g as 

IC:g^ g(y) dy. 

J — OO 

It is easy to see that d x JC = I and if lim^^-oo g{x) — 0, then ICd x = I. The problem here is that 
JCT Qca M s <£ L 2 (R). More precisely, while /CC c r a = ~Q ca G i 2 (R) we have that in general JCQ a L 2 (R). In 
fact, if f{u) — u p , then 

and therefore 

Since Q c is positive, /C£™ a is not an L 2 (R) function if p ^ 3. We remark that if p = 3, then there are no 
problems. This case is in the special class of nonlinearities considered in |2l)j . 

Our remedy to the above problem is to "regularize" the symplectic operator JC. Let Pq be the L 2 (M) 
orthogonal projection onto the subspace spanned by the translation vector and let Pq be its orthogonal 
complement. Then we define the anisotropic regularization /Cq q of IC as 

JC Qa :=KP Q +d- x P Q , 

where d' 1 :— (d x + a) -1 is the regularization of JC. We do not regularize in the direction of C*a since JC is 
well behaved on this vector. 
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For d a 1 to exist, the parameter a must lie in the resolvent set p{d x ) = C \ £R, and in such a case d a 1 
acts explicitly as 



on all L 2 (R) functions g. The lemma below, proven in Appendix iBl collects some properties of d a 1 , which 
will be used in the course of proving the main result. 

Lemma 3. Let 0, ip S L 1 (]R) n L 2 (R) and ael + . Then we have 



1. The operator d a l commutes with d x and spatial translation S a : f(x) — > f(x — a); that is, d x d. 
d^dx = I - ad^ 1 and Sad^ 1 = 9~ x 5 a . 

2. lid- 1 



-l 



| £ « < II^IIli 



3. There is a constant C such that ||^ a 1 ^|| i 2 — Get 2 \\ 
I Ifx<j) £ L l (R), then \\xd~ 1 <j>\\ L1 < C (art ||^|| 
5. If x4>, xip € L 1 (R) 7 i/ien | (0, c?" 1 ?/') — (0, A^) | < a (||<^||j,i H^IIl 1 + II^IIli II V'lliO an ^> in particular, 



II 1 - 



•a 2 so 



)0?X 



<2a||0|| £1 ||a;0|| £l 



6. Ifx^,x 2 ^> S L 1 (M) > f/ie/i 



dx + 01. 



For a small, the above lemma implies that the properties of K. and d a 1 are similar. Thus, we require in 
(fTTfl that 

^±JC Qa T Qca M s . (12) 

The existence and uniqueness of parameters a and c such that £ = u — Q ca satisfies <|12[) follows from the 
next lemma concerning a restriction of /Cq q and the implicit function theorem. 

The restriction Kq of K,Q a to the tangent space TQ ca M s is defined by the equation KqoVq = 
"PqICqo^Pq, where Vq is the orthogonal projection onto Tq cci M s . In the natural basis {Cta>Cca} °f the 
tangent space Tq oci M s , the matrix representation of Kq u is JV _1 f2Q Q , where 



._ | IICcoIIl 2 



IIC™JIl 2 



and 



Q ._ | (Ceo' K-QaCca) (Cca> K-QaCca) \ /io\ 

Notice that the matrix J7q q depends on the base point Q C a, and hence on the parameters a and c. 
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Lemma 4. If S'(c) > on a compact set I C and a << [S' \ :— inf/ 6', then Qqu is invertible for all 
c £ / and a £ R, and 

°£ = ^( KHoC^) 2 ^0+0° (u) 

g ^'(c) 2 V -s'(c) o / L*'J 

Hence, ||fi^|| = /or L<5'J smaZ/. 

Proo/. We use the relations C*a = "d x Q ca , lC Qa C c r a = ^QaC" a = 9 a 1 C? a i anti-self-adjointness of <9 X , 
and /C9 X = / to simplify the matrix £lQ a into 



fir 



-(CQca) 



Next, using statements ^ an d El of the previous Lemma, we separate the leading order part of flQ a from the 
higher order parts, and use the relation S'(c) = (Q ca , Q a ) to obtain that 

( ~ S ' (C) 2 W \ 

where | i? | < 2a sup/ HC^alli, 1 II^CSJIli- With Q ca and C™ a exponentially decaying, the estimate |(Q cq , c?" 1 ^) | < 

||Q ca || L i IKcnlLi is c l ear from tne properties of Thus, if a < \ \5'\ (sup 7 ||Q ca || L i ||C£zllz,i)~ i tncn tnc 
determinant 

det n Qa = 8\cf + aS'(c) (Q ca , d^Q a ) > ~ [S'\ 2 , 

and hence it is nonzero for all c £ I and a el. The coadjoint formula and the above estimate give 114fl . 
The estimate of ||^q^,|| follows from 114|) and the assumption that £™ a dx ^ 0. □ 

Given e > 0, define the tubular neighbourhood U e :— {u £ L 2 (R) | inf( c , a ) e / X R \\ u ~ QcaW^ < e} of the 
solitary wave manifold M s in L 2 (R). 

Proposition 5. Let I C M+ be a compact interval such that c i— ► Q ca is C' 1 (J). Then, if a << [S'\, there 
exists a positive number e = e(I) = 0a 3 \_5'\ 4 arid unique C 1 functions a : U e — ► M+ and c : U e — ► J, 
dependent on a and I , such that 

(Qc{u)a(u) ~ u ^QaC c \u)a{u)) = a « d (<2c(u)a(«) ~ u ^QaC(u)a(u)) = 

/or a/? u £ U e . Moreover, there is a positive real number C = C(I) such that 

\\u~ Qc(u)a(u)\\ H i < Ca~^ inf^ ||u - Q oa || H i (15) 

/or aZ/ue ^nF 1 ^). 
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Proof. Let fi := (a c) T and define F : L 2 (R) xR + xJ^B 2 by 

The proposition is equivalent to solving F(u,g(u)) = for a C 1 function <?. Observe that F is C 1 and 
F(Q ca , /i)=0. To apply the implicit function theorem it suffices to check that d p F(Q cai fi) is invertible. 
Then the there exists an open ball B e (Q ca ) of radius e with centre Q ca , and a unique function g a Q ca ■ 
BeiQca) -> K + x 7, such that F(u, g aQaa (u)) = for all w G B E (Q ca ). Since d li F{Q ca ,n) = VL Qa , the 
invertibility of d^F(Q ca , fi) follows from Lemma 0] provided a is small enough. The radius of the balls 
B e (Q ca ) depend on the parameters c, a and a. To obtain an estimate of the radius, and to show that we 
can take e independent of the parameters c and a, we give a proof of the existence of the functions g a Q aa 
using the contraction mapping principle (just as in the proof of the implicit function theorem). 

Expand F{u,fi) to linear order in fi around /io = (ac) T : 

F(u, fi) = F(u, no) + d p F(u, fi )(fi - /i ) + R(u, fi), (16) 

where R(u, fi) = ^d 2 F(u, (1 — A)/i + A/i)(/i — /io) 2 for some A € [0, 1]. The operator d p F{u, fi ) is computed 
to be d^Fiu, fi ) = Oq q + A, where 

a ( {Qca - U, d a Qca) (Qca ~ U, d c Q ca ) \ 

■ \ (Qca - U, dad^Qa) (Qca ~ U, d^Qa) • / 

If it G B £ (Q ca ), where e remains to be chosen, then the properties of d~ x imply that ||A|| < Ca~^e. Thus, 
since flQ a is invertible, if e « (C sup/ H^qqJI) -1 ^ 3 , then d^F(u,fJLo) is invertible and || /io)] -1 1 < 

C sup/ || Oq^, (I . Hence, given u, F(u, fi) — has a solution fi if and only if 

fi = H(fi) := /i - [dfj,F(u, /io)] -1 (F(u, fi ) + R(u, fi)) 

has a solution /i. The latter is equivalent to the function H having a fixed point. This is guaranteed by the 
contraction mapping principle if H is a strict contraction from some ball B p (fi ) to B p (fi Q ). 

Say fi G -B p (/io), where p remains to be chosen, and consider the bound 

||ff(M)-Mo|| <Csup||n^||||F(u,Mo) + J R(n,M)ll- 

After subtracting F(Q ca , (Iq) — from F(u,/io) and using the mean value theorem, the above becomes 

||F(/i)-Mo|| < Csup||0 Q i|| (\\d u F((l-X 1 )Q ca + X 1 u,fi )\\ \\u-Qca\\ 
I 

+i ||d 2 F( M , (1 - X 2 )fi + X 2 fi)\\ \\fi - /io|| 2 ) 
for some Ai, A 2 € [0, 1]. Again, using the properties of d~ x we find that 

\\d u F({l - Ax)Q ca + Am, Mo)|| < Ca-i 

\\dlF{u, (1 - A 2 )mo + A 2 /i|| < C(l + a - ^) (17) 
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for all and u £ B £ (Q ca ). Thus, if p < 1, then ||-ff(^) - jUo|| < Csup z ||0g^,|| (aT5£ + p 2 J. 

Taking e << (C sup/ ||S1q*| )~ 1 a^p and p << (C supj- ||^Qa||) _1 implies H maps B p (p ) into B p (p ). 
Let pi,p2 S Bp(fio) and consider the bound 

||#( M2 ) - < Csup||n^|| ||i?(u ;Al2 ) - i?(u, m )|| 

or, using the mean value theorem, the derivative of I jlfijl with respect to /i, and the mean value theorem 
again, 

\\H(jm 3 ) - H(fxi)\\ < Csup ||fi^|| \\dlF(u, (1 - A 2 )[(l - AO/xi + A lA *a] + A 2 /i )|| 

x ||(1 - Ai)/ii + Ai/i 2 - Mo|| \\fi2 - pi\\ 
for some Ai, A 2 £ [0, 1]. Using l(T7|l and ||(1 — \\)p\ + X1P2 — Moll < P then gives 

||ff(/x 2 ) - fT(A*i)|| < Csup||0^||(l + a~h)p\\p 2 - p x \\. 
1 

Thus, with the above choices of e and p, H is a strict contraction. We conclude that the radii of the balls 
B e (Q ca ) can be taken independent of a, c (but dependent on /) and e = &ah [5'(c)J 4 . 

The above argument shows that there exists balls {B £ (Q ca ) \ a £ R+,c £ /} with radius e dependent 
only on the parameter a and the compact set /. Notice that U e = (J{ B £ (Q ca ) | a £ R + , c € I}. Pasting the 
C 1 functions g a Q aa together into a C 1 function g a j : f/ E — > K + x / gives the required C 1 functions a(it) and 
c(it). Uniqueness follows from the uniqueness of each of the functions g a Q ca - 

Let u £ C/ E ) c£ / and o£l, and consider the equation 

u Q c(u)a(u) — u Qca + Q ca Qc(u)a(u)- 

Clearly, inequality (fT3|l will follow if ||Q ca — Q c (u)a(u)\\ H i < C\\u— Q ca \\ H i for some positive constant C. 
Since the derivatives d c Q ca and d a Q C a are uniformly bounded in iJ 1 (R) over / xl, the mean value theorem 
gives that ||<5 ca — Q c (u)a(u) \\ H i < C||( ac ) T — ( a ( M ) C ( U )) T ||: where the constant C does not depend on c, 
a, or a. The relations g a i(Q ca ) = {ac) T and g a i(u) = (a(u)c(u)) T then imply \\Q ca ~ Qc(u)a(u)\\ H i < 
C||ffa/(Qca) ~3a/(w)||. Again, we appeal to the mean value theorem and obtain (|15fl . using the properties of 
9" 1 and that d u g a i = —d^F~ 1 d u F is uniformly bounded in the parameters c £ I, a £ K + and w g f/ E . □ 

4 Spectral Properties of the Hessian d 2 A ca 

The Hessian d 2 A ca at Q ca in the L 2 (K) pairing is computed to be the unbounded operator 

C Q :=-d 2 x + c- f{Q ca ), (18) 
defined on L 2 (R) with domain H 2 {M). We extend this operator to the corresponding complex spaces. 
Proposition 6. The self-adjoint operator Cq has the following properties. 
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1. C Q & = and C Q Q a = -Q ca . 

2. All eigenvalues of Cq are simple, and Null Cq = Span {C*al- 

3. Cq has exactly one negative eigenvalue. 
4- The essential spectrum is [c, oo) C R+. 

5. Cq has a finite number of eigenvalues in (— oo, c). 

Proof. Recall that the vectors C*a := —d x Qca and Q a :— d c Q C a are in the Sobolev space H 2 (R). Thus, 
relations CqC^ — and Cq^ — —Qca make sense, and are obtained by differentiating A' ca (Q ca ) = with 
respect to a and c. The first relation above proves that is a null vector. 

Say £, n 6 ff 2 (R) are linearly independent eigenvectors of £q with the same eigenvalue. Then, since Cq 
is a second order linear differential operator without a first order derivative, the Wronskian 

is a non-zero constant. With rj and £ both in iJ 2 (R) however, the limit lmx^oo W(rj, £) is zero. This 
contradicts the non vanishing of the Wronskian, and hence all eigenvalues of Cq are simple and, in particular, 
Null C Q = Span 

Next we prove that the operator Cq has exactly one negative eigenvalue using Sturm-Liouville theory on 
an infinite interval. Recall that the solitary wave Q C a(x) is a differentiable function, symmetric about x = a 
and monotonically decreasing if x > a. This implies that the null vector or equivalently, the derivative 
of Qca with respect to x, has exactly one root at x = a. Therefore, by Sturm-Liouville theory, zero is the 
second eigenvalue and there is exactly one negative eigenvalue. 

We use standard methods to compute the essential spectrum. Since the function f'{Q C a{x)) is continuous 
and decays to zero at infinity, the bottom of the essential spectrum begins at lim :E _» 00 (c — f'{Q C a{x))) = c 
and extends to infinity: a ess (CQ) = [c, oo). Furthermore, the bottom of the essential spectrum is not an 
accumulation point of the discrete spectrum since f'(Q C a(x)) decays faster than x~ 2 at infinity. Hence, there 
is at most a finite number of eigenvalues in the interval (— oo, c). For details see . □ 



5 Anisotropic Coercivity of the Hessian Cq on {KLQ a TQ ca M s ) 

In this section we prove strict positivity of the Hessian Cq on the orthogonal complement of the 2-dimensional 
space K,Q a TQ aa M s = Span {Q ca , d^ 1 ^}- This result is a crucial ingredient in the proof of the bound on 
the fluctuation £. 

Proposition 7. Assume S'(c) > on the compact set I C To- The following statements hold if a > is 
small enough and ^_L/CQ Q Tg ca M s . 

1. There are positive real numbers C\ and C2, independent of a, and a function q{ol) satisfying C\u < 
g(a) < C2CX such that (Cq^,^) > g(a) \\£,\\ H i for all c G I and agl. 

2. The infimum inf{(£g£, £) \^l.K,Q a TQ ca M s and ||£||_l2 = 1} * s attained and the unique minimizer r\ is 
of the form r\ = jCca + V±! where \\n±\\ H i = Q)ak and 7 = 01. 
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3. Let rj be as above and let £ g := £ — (£, rj) n. There exists a positive real number C3 independent of a, 
such that (£q£o,£o) > C3 ||£ 9 ||jji (notice that a enters £ s through the minimizer rj). 

4- The Hessian Cq is anisotropically coercive on (K.Q a TQ ca M s )~ L ; that is, 

(C Q ^0>C3U a \\ 2 ^+CaU b \\ 2 H1 . 



Proof. Define the set X := {£ £ iJ 1 (R) | ^±ICQ a TQ ca M s and ||£|| i2 = 1}- Our first step is to prove an upper 
bound on infx (£q£,£). We do this by computing (£q£, £) for the test function 

£ := AiC*a + ^2^~ Qa + A3<3ca, 

where Ai, A2, and A3 are chosen to satisfy (£, Q ca ) = 0, (£, c^Coo) = ^' ano - II£Hl 2 = These conditions 
imply £ € X, and, after substituting £ with its definition, have the form 



A 2 (Q C a,9- 1 C) + A 3 ||Q ca | 



2 

L 2 
l/-n 



Al (Cca, 1 (ca) + A 2 ||<9 Q 1 Cca|| L 2 + A3 (Qca, d a 1 Q a ) 



and 



A? ||Cca|| L 2 + ^2 ||9 a 1 Cca 11^2 + A3 ||Qco|Il2 



+ 2AiA 2 (C d-'Q,) + 2A 2 A 3 (Q ca , CO = 1. 



(19) 
(20) 



(21) 



Equation 1)21) [) can be solved for Ai when (Cca ^ a 1 Cca) * s n °t zero. A straightforward computation using 
antisymmetry of d x and statements ^and [5] of Lemma [21 gives that 



(C,d- 1 C a ) = 6'(c) + 0c 



(22) 



Thus, if a is small enough, then <9 a 1 C™a) * s non-zero. We substitute for A x in lf?T|l using lt2U|) and then 
use ifH?)) to substitute for A 3 to obtain 



IC 



ti-||2 
call L 2 



\Sca,°'c< S> CQ / 



a Si 



Cca) 
11*3 ca|| £2 



l<9 C 



a Sca||£2 



IIQcll^ 



A^ = l. 



(23) 



This equation, the relations ||5- 1 C" a || i2 = £ C"a + 01, (Qca, CCS*) = 01, and © imply that 

A 2 = 0a. Equation (|19fl then implies A3 = 0a. Evaluating the quadratic form (Cq-, •) at the test function 
£ and bounding with Holder's inequality implies 

< C(c)a. 

In the last inequality we have used the bounds on Ai and A 2 , and the above estimate of c^Coollfl-i- The 
constant C(c) does not depend on the parameter a since H 1 (M) and L°°(R) norms are translation invariant. 

To prove the first part of the proposition, we first prove that infx (CqCO > 0, or equivalently, that 
inf x nff 2 (R) (Cq£, £) > 0. By the max-min principle, infxn_f/ 2 (R) (Cq£, £) is attained or is equal to the bottom 
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of the essential spectrum. We take a small enough so that the above upper bound is below the essential 
spectrum. Let 77 be the minimizcr. 

We claim the set of vectors {C*ai Cccn v} is a linearly independent set. If they were dependent, then, since 
£*a and (ca are orthogonal, there are non-zero constants 71 and 72 such that 77 — 7iC*a + liCca- Projecting 
this equation onto Q ca and d' 1 ^ gives the equations j 2 S'(c) = and 71 d^Qa) +72 (d^Qai C*a) = °- 
Thus, if a is sufficiently small, then both constants are zero (we have used l|22|l). This is a contradiction 
since the zero function does not lie in X . 

The above argument proves that Span {C'a; (cm v} is three dimensional. The min-max principle states 
that if 

E 3 := inf sup (£q£,Q 

VCK 2 (R),dimV=3 SeV) || ? || t2=1 

< max (£n£,£) 

£eSpan{C*S,C a ,77},|ISIIi,2=l 

is below the essential spectrum, then it is the third eigenvalue counting multiplicity. Let £ = 71 77+72(^2 +73 C*a 
be the maximizer in the second line above. By Proposition there are exactly two non-positive eigenvalues. 
Hence 

0<E 3 < (C Q £, £) = 7l 2 (C Q r h 77) - jIS'(c). 

Thus, since 5'(c) > 0, we must have (CQr],r]) > 0. The function <r(c, a) :— {CQT],ri) is continuous with 
respect to c since both Q ca and c*" 1 ^ are continuous as mappings taking c to elements of Taking 
the infimum of <r(c, a) over / implies (£q£, £) > p( a ) \\£,\\l 2 f° r au C-L^Qa^Q C a-^«) wnere p( a ) := mi / c( c > a )- 

To complete the proof of the first statement (modulo the lower bound on g(a)), we improve the above 
lower bound to one involving ff 1 (R) norms. If we define K :— sup 7 ||c — /'(Q CQ )|| L=0 , then (£q£,£) > 
\\dx£\\ 2 L 2 — K \\£\\ L 2 for all £ € iJ 1 (R). Adding a factor -^ij of this bound to the above bound gives the 
required result with 

p(a) 
p(a) + A + 1 

Notice that the upper bound a(c,a) < C(c)a derived above gives, after maximizing constants over c G /, 
the uniform upper bound on g. 

As already shown, the minimizer 77 of infx (£q£,£) exists for a small enough. We prove the properties 
of 77 by manipulating its Euler-Lagrange equation 

£ Q 77 = 0r, + (3iQ ca + (24) 

where (3, pi, and /3 2 are Lagrange multipliers for the constraints ||?7|| L 2 = 1, (77, Q ca ) = 0, and (77, d" 1 ^) = 0. 
The inner product of this equation with 77 shows that f3 = a. Take a small enough so that (3 — a is not 
an eigenvalue of Cq. Then the minimizer is unique. Indeed, the difference ( between two minimizers is a 
solution to Cq( = P(. Since (3 is not an eigenvalue, ( = is the only solution to this equation. 

We now decompose 77 orthogonally as 77 = 7^ + V± ! substitute this decomposition into Q24J1 , and use 
^-qCcq = to obtain the equation 

(Cq - a) 77x = lo£ a + piQca + /WC (25) 
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for rjx- To solve for rjx, we first project this equation by Pq, where Pq = 1 — Pq and Pq is the orthogonal 
projection onto the nullspace of the operator Cq: 

[Lq - a] V± = PxQca + faPqd^Q,, 

where Lq is restriction of Cq onto the orthogonal complement of the null space of Cq. The spectrum of Cq 
has essential spectrum [c, oo) and a finite number of eigenvalues in (— oo,c). Moreover, Cq is independent 
of a, and therefore, for a small enough, we conclude that the interval [0, a] is disjoint from the spectrum 
of Lq, with the distance to the spectrum of Lq bounded below by a positive number C independent of a. 
Hence, we can solve (12511 for rjx to obtain 

r,x = (L Q - ay 1 (faQ ca + frpQ&^Q) . (26) 




To prove the L 2 (R) the estimate of rjx, we require estimates on the Lagrange multipliers fa and fa. We 
take the inner product of i|25[l with and use (|22|) to obtain that 

fa[5'{c) + 0a] = - 1 a\\£ l f L2 . 

Thus, fa — 0a since the constraint H^H^a = 1 implies 7 = 01. Similarly, since r\ is orthogonal to Q ca , the 
inner product of l|24|) with £™ a and statement [S] of Lemma [3] gives the relation 

= -<r<0,Ga)-JhS'(c). 

The estimate (3\ = 0a is immediate using the estimate of P2 and the assumption <5'(c) > 0. We substitute 
the estimates of fli and fa, and the estimate ||9~ 1 ('™ a || i2 = 0(a _ 2) into (|26|l . and use the above fact that 
a is at least a distance C away from the spectrum of Lq to conclude that 

\\nx\\v = 0a-ia. (27) 

Replacing a with its upper bound gives the third statement of the proposition. 

We now prove a lower bound on the infimum a = infx (Cq^,^). We again need to take a small 
enough so that a minimizer exists. The orthogonal decomposition rj = "/Cca + V-L °f the minimizer implies 
a = (CQr/±,r/±). Substituting for Cqt/± using (|25|l gives 

v = v\\vAli+fc{rU,d- 1 £ a ), 

where we have used that r/j_) = and, since the minimizer r/ is orthogonal to Q ca , (rj±, Q ca ) = 0. Thus, 
IP7||. Bn = 0a, and ||9- 1 Cca|| i 2 = 0" _1 imply 

a 3 a 2 
a <d— + C 2 — 

a a 

or, since a is positive, C\a 2 + C 2 a — a > 0, where the constants C\ and C 2 depend continuously on c. The 
positive root of the quadratic is a lower bound on a. After rationalizing, we obtain 

2a 

° > , , > Kia, 

C 2 + VC% + 4 ^ia 
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for some constant K\. Minimizing the constant over / completes the proof of the lower bound. 

Our proof of statement three in the proposition requires that infy- is positive, where Y := {£ € 

i/ 1 (IR) | C-LQc C*a an< i II£IIl 2 = ■*-}■ The argument is similar to the proof of inf x (£q£, £). By the min-max 
principle either infy n H 2 (R) = infcr ess (£Q) or the minimizer is attained. There is nothing to prove in the 
former case since inf <r ess (£Q) = c > 0. Thus, we assume 77 is a minimizer. As above, the set {77, C™ a } 
is a linearly independent set due to the assumption S'(c) > 0, and the min-max principle implies the third 
eigenvalue E 3 satisfies 

0<E 3 <7i{£ Q r),v) 

for some constant 73. Thus, we must have (£q£,£) > K 3 for all £ € Y, where K3 is a positive constant 
independent of a. As with the infimum over the set X, this inequality can be improved to the H (K) estimate 
(C Q t,Q>C 3 \\tf B1 forall£lQ CO) £. 

We now decompose £ g orthogonally as £ ff = + ip. Since £ 9 is orthogonal to -q and r\ = 7^ + ?7-Lj 

= IKcall^T" 1 (tg,V-Tlx) = 0\\VX\\V MbWi? = UAv 

Substituting this bound into ||^|| 2 ffl > ||^||^ - /3 2 O^i gives that > ||£ 9 ||^ (1 - a K%\f m ). 

Thus, if a < h WClaWn 1 ' then H^H^i > \ Substituting this into the inequality 

(£Qt g ,Q = {£Qi>,i>)>c 3 U g f H i 

(which follows from the fact that ipl. Q ca , C*a) completes the proof. 

To prove the last statement we define £{, := £ — £ g . Since the vectors £ and rj are both symplectically 
orthogonal to the tangent space, so is Thus, using the above inequalities, 

(CqCO > C 3 \\t g \\ 2 H1 +Ca\\t b f H1 +2{£,r } ) (C Q r,,Q ■ 

The cross term (Cqij, £ s ) is zero; indeed, substitute for Cqt] using equation (|24|l and use (£ 9 , 77) = to obtain 

(C Q r,, Q = ft (Q m , Q + 2 (9"^, £„) . 

This expression, however, is zero since both £ and rj are orthogonal to Q ca and d~ l C,^ a , completing the 
proof. □ 

6 Evolution Equations for the Fluctuation £ and the Parameters 
a and c 

In Section we proved that if u remains close enough to the solitary wave manifold M s , then we can write 
a solution u to Q uniquely as a sum of a modulated solitary wave Q ca and a fluctuation £ satisfying 
orthogonality condition (|12fl . Thus, as u evolves according to the initial value problem QJ, the parameters 
a(t) and c(t) trace out a path in R 2 . The goal of this section is to derive the dynamical equations for the 
parameters a and c, and the fluctuation £. We obtain such equations by substituting the decomposition 
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u = Qca + £ into Q and then projecting the resulting equation onto appropriate directions, with the intent 
of using the orthogonality condition on £. 

From now on, u is the solution of with initial condition u$ satisfying infQ ca eM« \\uo — QcaWjj 1 e ' 
and Tq = Tq(uo) is the maximal time such that u(t) € U e for < t < Tq. Then for < t < Tq, u can be 
decomposed as in <|11|) and \1'2\, 

The majority of the work involves estimating the higher order terms of the resulting equation for the 
modulation parameters. It turns out that a naive attempt at bounding £ directly with the Lyapunov method 
does not give good results. As will be seen later, the component of £ in the direction of (* a is particularly 
problematic. On the other hand, C, l ca = —d x Q ca is the derivative of a function and the null vector of Lq. 
This can be used to improve the bound. Thus, in order to obtain better estimates on £, we orthogonally 
decompose the fluctuation as 

£ = & + 

where = (£, 77)77. Recall that 77 is approximately and is given in Section [5j 77 = jC^a + r]±, with 
\\r]\\ L 2 — 1, 7 = 01, and ||?7±|| £ 2 = 0a 5 - We use the above decomposition to prove the following proposition 
regarding the dynamical equations for a and c. 

Proposition 8. Assume S'(c) > on the compact set I C K+. Say u — Q ca + £, is a solution to l|T]l. where £ 
satisfies (|12f) flTid £ = + £ g as above. If oT^ \\£g\\ H i + ||£||#i is small enough and a, e x < 1, then, provided 
eel, 

0) - (-r>) + ^(-i(nja*)') + ^o. (»> 

iu/iere |Z(a, c,£)| < C (ae a e x + e a e^ + (a? + a~^e a e ;E ) Ug\\ H i + (a + e a e x ) \\£\\ H i + U\\%i) , for some pos- 
itive constant C = C(J). 

Proof. Recall that the solitary wave Q ca is an extremal of the functional A CQ . To use this fact we rearrange 
definition JSJ of A ca to write the Hamiltonian H D as 

1 f°° 

H b {u) = A ca (u) - cP(u) + - bu 2 (x) dx, 

* J -00 

where for notational simplicity we have suppressed the space and time dependency of b. Substituting Q ca + £ 
for u in and using the above expression for Hi, gives the equation 

aCZ + cCa + £ - d x A' ca (Q ca + 0- cd x [Q ca + + d x [(Q ca + £)b], 

where dots indicate time differentiation. Taylor expanding A' ca (Q ca + £) to linear order in £ and using that 
Q ca is an extremal of A ca gives 

£ = d x [(C Q + * + 6(a) - c)£] + d x N'(0 - [a - c+ 6(a)]C - 

+ 6'(a)a x [(a; - a)Q ca ] + d x [S 2 bQ ca }. (29) 

We have used the relation = —d x Q ca , definition (|18|l of £q, the definitions 

* := b(x) - b{a) 
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and 

S 2 b := b(x) - b{a) - b'(a)(x - a), 

and definition l|5U|) of N'(£) given in Appendix El to write the above equation in a convenient form. Define 
the vectors £1 := C*a an d C2 := C"<r Projecting l)29|) onto K.Q a Q, for i = 1 and 2, and using the anti-sclf- 
adjointness of 9 X gives the two equations 

[d - c + 6(a)] [<C /C Qq C,) + & S s /C Qa Ci)] + c <C, K Qa Q + (£, K Qa Ci) - a (£, d x K,Q a C,i) = 

- b'(a) ((x - a)Q ca , d x K, Qa (i) - (S 2 bQ ca , d x JC Qa C,i) 

- (£ Q £, SWCgaCi) ~ ^/C Qa Ci) - &/Cq q G) • (30) 

We can replace the term containing £ since the time derivative of the orthogonality condition (£, /Cq q £j) = 
implies (i,K. Qa Ci) = a (£,d x K,Q a (i) - c{^dJC Qa Ci)- Note that we have used the relation d a Q = -d x Q. 
Thus, equations i|3U|) in matrix form are 

(I + B)n Qa ( h - C t bia) )=X + Y, (31) 



c 



where 



' : 1 (£ Q £ + *e + 7V'(£),C-a9- 1 C) 

and 

(£, Cca) <C,C 



* —™ ( ) + >'(°> ( «<(*-«)l,^> ) - ( 



fl := 



(CO -(C^C*) 



We have explicitly computed ((x — a)Q ca ,Q) and used the relations d x K.Q a Cca = Cca an d ^x^QaCca = 
C" Q ^ ac?" 1 ^ to simplify the above expressions. 

We now estimate the error terms and solve for d and c. The assumptions we made on the potential 
imply that 

|*| < e a e x (x - a) and \5 2 b\ < e a e 2 x (x - a) 2 . (32) 
Thus, by Holder's inequality and the L°°(R) estimate \\d~ X Q || < ||Ci||j,i 5 

||X|| - -b'(a)6(c) f J J + ®^e a e x + e a e 2 x 
= &e a e x . 

In the last equality we have used a < 1 to bound ae a £ x by e a £ x - We now estimate the vector Y using 
the properties of Cq given in Appendix {7\ Indeed, the generalized nullspace relations CqC^ = and 
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CqCcu — ®xCca i m ply £i) = 0, where the orthogonality condition (£,^QaCca) = is used when i = 2. 

Although the same is not true of the inner product a (Cq£, c*" 1 ^), we still use the relation £qCco = to 
obtain a bound. If we use rj = 7^ + r]± , then the orthogonal decomposition + £ g of £ becomes 

where 7 = 01, ||r7j_|| L2 = 0a 5 and ||r/|| i2 = 1- Thus, 

a (C Q t d-'Qa) = « (Cq^B- 1 ^) + a (£, rj) (£ Q 7?±, c^C) 
= ®^\\i g \\ Hl +a\\i\\ H1 . 

The terms containing the potential, but not the term d~ l C£ a , are easily estimated using the bound on <55 and 
exponential decay of and Q a . The resulting estimates are 

= M\\ La and (*£,C) = @e a e x ||£|| L2 . 

To obtain a bound of a <9~ 1 CSj/> we decompose £ as above to obtain 

Then the estimates 1 1 1 CSi 1 1 = 01 and ||x9~ 1 Cca|| L 2 = 0a~^ imply 

Lastly, Holder's inequality and Lemma El imply (iV'(£),£i) = 0||£||^i and a (N'(£), d~ l Q a ) = 0a? U\\ 2 m . 
Adding the above estimates gives 

11*1 = + oT^cex) II^H^ + (a + e^) ||£|| ffl + . 

The inner products (£, C*a) an d (£>C"a) are clearly of order We estimate the remaining entry 

(£, dcd^Qa) = (Cj d^dlQca) of i? by the same technique as above: we replace £ by 7 (£, 77) C*o+Cg + (^ »?) »7-L 
and estimate the resulting expression to obtain 

P*||=0cr* + UWm- (33) 

We take ||S|| smaller than one so that 7 + B is invertible and ||(7 + 1 1 1 = 01. Acting on (|31() by 
(I + B)- 1 = I - B(I + B)- 1 and then Q Qa gives that 

a — c + 6(a) 



c 



= r> Qa x - n Qa B(i + b)- 1 x + Q Qa (I + B)~ 1 Y. 



We use the leading order expressions for £tQ a and A, and the bounds on ||A||, ||F||, \\B\\ and ||(7 + B) x | 
to obtain the estimate 



V c / S'{cf \ 5 > (c) 



6>(c) 

+0ae a e x + e a e 2 x + ( Q i + oT^ex) ||f fl || H , + (a + e a e x ) ||f + . 

In the order notation used above, the implicit constants are continuous with respect to the parameter 
c and independent of the parameter a. Maximizing these constants over the compact set I completes the 
proof. □ 
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7 The Lyapunov Function 



In the last section we derived dynamical equations for the modulation parameters. These equations contain 
the H 1 (R) norm of the fluctuation. In this section we begin to prove a bound on £. Recall that the latter 
bound is needed to ensure that u remains close to the manifold of solitary waves M a for long time. 

We employ a Lyapunov argument with Lyapunov function 

M c (t) := A ca (Q ca + - A CQ (Q ca ) + b'(a) ((x - a)Q ca , £} . (34) 

Remark: if f(u) — w 3 , the last term in the Lyapunov functional is not needed; however, apart from compu- 
tational complexity, there is no disadvantage in using the above function for this special case as well. 

Lemma 9. Assume 8' (c) > on the set I . Let u = Q ca + ^ be as in Proposition^ Let a, e a , e x , ||£|| H i < 1. 
Then, if oT^ + ||£||#i * s sma ^ enough and c G I, there is a constant C — C(I) such that 

^M c (i)<C[( £a e K + ||£|| 2 ffl )S], (35) 

where E := e a e x + (a? + a-h a e x ) U g \\ H1 + (a + e x + e t ) \\£\\ H i + U\\ 2 H i- 



Proof. Suppressing explicit dependence on x and t, we have by definition 

1 



A CQ (it) := H b (u) 



u bdx + cP(u). 



Thus, relations JIJl, © and JO) imply that the time derivative of A CQ along the solution u is 



dt 



Ka{u) 



l 2 +b' 



l C u 2 -uf(u) + ^(d x u) 2 + F(u) 



+ b" ud x u dx. 



Substituting Q ca + £ for u, manipulating the result using antisymmetry of d x , and collecting appropriate 
terms into b'(a) (£ Q £, d x ((x - a)Q ca )), (N' (0,d x [8b(Q ca + £)}) , and (A' ca (Q ca ), d x (Sb(Q ca + £))) gives the 
relation 

^[A ca (Q ca + f) - A ca (Q CQ )] =b'(a) (C Q ^d x ((x - a)Q ca )) +c(Q ca ,Z) + {C Q ^d x (5 2 bQ ca )) + ||£||* a 

+ 4 ( b 'Z> + | ( & '^, d x - (f'(Q ca )t d x ($>0) 

+ (jV'(£), d x [5b(Q ca + £)]) + (&"£, d x + (A' ca (Q ca ),d x [S>(Q ca + £)]) • 

The last term is zero because A' ca (Q ca ) = 0. The inner product (£, Q c ) is equal to —a^,d~ 1 Q ca ) = 
H x since -d^Cca = Q ca - ad a 1 Q ca and £,±Q ca - We use Lemma El assumptions |JTDJ on the 
potential, estimates and 



\Sb'\ < e a e 2 x x 
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to estimate the size of the time derivative. We also use that Q ca , d x Q ca , d x Q ca and f'(Q ca ) are exponentially 
decaying. When e x < 1, higher order terms like (6"£, d x £} are bounded above by lower order terms like (6'£, £) . 
Similarly, if < 1, then e a e x ||£||#i < e a £x H^llffi- This procedure gives the estimate 

J t [Ka{Qca +0 - AcaiQca)} =&'(<l) (£, C Q d x ((x - a)Q ca )) + (iV (£) , <56cU) 

+ 0|c|||€||ff 1 +e«4lklliP+c«c*IICII ^i. 

We compute 



^ (V(Q ca + - F(Q ca ) - f(Q ca )£ - \f'{Q ca )e^j b' dx, 



and use the second estimate and the proof of the third estimate of Lemma IT21 to obtain (N'(£), Sbd x £) = 
0e a e x U\\ H i- Thus, since e a e x \\£\\ H i < e a e x ||£|| H i when \\£\\ H i < 1, we have 

j t [Aca(Q ca + - A CQ (Q ca )] = 6' (a) (£, /^((a; - a)Q ca )) 

+ 0|c| lieilffi + e a £ IICII^ + e Q e, ■ (36) 

When /(u) = u 3 , (£,£Qd x ((x — a)Q ca )) = since £™ Q = d x [(x — a)Q ca \. In this special case the above 
estimate is sufficient for our purposes, but in general, we need to use the corrected Lyapunov functional. 
When £ e C(R, H 1 ^)) n C 1 (R, #~ 2 (1R)), 6' (a) (£, (as - a)Q ca ) is continuously differentiable with respect to 
time; 

j t [b'(a) (£, (a - a)Q ca )} =d t b' (£, (x - a)Q ca ) + &'(a) (£, (x - a)Q ca ) + c&'(a) (£, (a - a)C> 

+ db'(a) <£, (x - a)C) + db"(a) (£, (x - a)Q ca ) , 
where (£, Q ca ) = has been used to simplify the derivative. Substituting for <9 t £ using l|29[) gives 

j t [b'(a) (£, (a: - a)Q ca )] = - 6' (a) (£, £ Q ^((x - a)Q ca )) - [a - c + b(a)]b' (a)± \\Q ca \\ 2 L2 + d t b' (£, (x - a)Q ca ) 

+ [d-c + b(a)]b'{a) {d x t, (x - a)Q ca ) + [a - c + 6(a)]6"(a) (£, (x - a)Q ca > 
+ c6'(a) (£, (x - a)C) - 6'(a) (£, - a )Q ca )) - b'(a) (N'(0, d x {(x - a)Q ca )) 

- b'(a) (5 2 bQ ca , d x ((x - a)Q ca )) + [c - 6(a)]6"(a) (£, (x - a)Q ca ) . 
We estimate using the same assumptions used to derive (|36H . If ||£||^i and e x are less than 1, then 

j t [b\a) (£, (x - a)Q ca )} = - b'(a) (£, £ Q ^((x - a)Q ca )) + 0|d - c + 6(a)|e Qe;E + |c|e ae;c 

+ 0e«4 + ((1 + £a)4 + e«et)c« ff i + ||£||^ ■ 

Adding the above expression to (|36(l gives an upper bound containing |c| and |d — c + 6(a) |. Replacing these 
quantities using the bound 



\c\ + \a-c+b(a)\ = 0e a e x + (a* +a h a e x j Ug\\ H i + {a + e a e x ) \\£\\ ffl + 



2 
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from Proposition [S] and bounding higher order terms by lower order terms gives <|35[) . To use the above 
bounds on |c| and |d — c + b(o>)\ we must assume WCgWm + is small enough so that Proposition [S] 

holds. □ 



8 Bound on the Fluctuation and Proof of Main Theorem 

We are now in a position to prove the bound on £. 

Proposition 10. Say u = Q ca + £ is a solution to Q], where £ satisfies Let e ai e x < 1 and < s < i. 

Then, if e a e x is small enough, there are constants G\, Ci, and C3, such that if the initial condition uq 
satisfies m{Q ca£Ms \\uq - Qca|| ff i << (e Q ex) 2s , then 

U(t)\\m <Ci(e a e x ) s and U 9 (t)\\ H1 < C 2 (e e x )*, 

for all times t < T\ :— C3 ((e a e x ) s + e t + e x ) 1 . 

Proof. We choose eo := miQ cae M B \\uq — QcaWa 1 smau enough so that ||£(0)|| H1 < Ca-2e a < e (see JT5J) 
is small enough to satisfy the conditions of Proposition [5] and Lemma Then, continuity of the solution 
u = Qca + £ m H (WL) with respect to time implies the conditions continue to be satisfied over a non-empty 
time interval [0, T]. We will obtain an estimate of ||£(t)||#i over a time interval [0, T] by deriving an equality 
for ||£(i)Hiji from upper and lower bounds on the Lyapunov functional. We suppress dependence on t for 
notational convenience. 

Define := sup[ T ] anc ^ l^l T : ~ su P[o t] 1^1- Integrating the time maximized upper bound 

in Lemma El gives 

M c (t) < \M C (0)\ + C{e a e x + )|3| T T 

for all t € [0, T] . A lower bound is obtained by expanding the A ca (Q ca + £) term in M c (t) to quadratic order 
and using A' ca (Q ca ) = to obtain 

M c (t) = {C Q H, + N(0 + b'(a) (£, (x - a)Q ca ) , 

where the nonlinear remainder N(£) is defined in Appendix IH1 Estimating ./V(£) with Lemma ll 21 and using 
anisotropic coercivity of the Hessian Cq (Proposition 0) gives 

M c {t) > C 3 U 9 \\m + Ca U b \\ 2 H1 - C(e a e x U\\ H1 + Uf m ). 
Together with the upper bound on M c (t), this bound implies 

U g f m , a ||6||^ < MgW* + <* MbWm < \ M c(0)\ + (t a e x + Uf T mTT + e a e x ||£|| T + Uf T . (37) 

Note that we have set non-essential constants to unity. Since the above inequalities hold for all t S [0, T], 
we can replace ||£b||jji with ||£h|j T and ||£ ff ||rxi with ||£ S || T . Multiplying the resulting inequality for ||£ ff || T by 
a and adding to the inequality for a W^bWrn gives 

a U\\t < \M C (0)\ + (e a e x + U\&)\S\ T T + e a e x U\\ T + U\\t ■ 
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Next, we take |^|tT = 0« and ||£|| T = 0a to obtain the bound 

U\\ T < a-i\M c (Q)\i + {e a e x )? + e^.a" 1 . 

The initial value of the Lyapunov functional M c (0) can be bounded by the iJ 1 (R) norm of the initial 
fluctuation ||£(0)|| ff i < Ca"2£ . Indeed, Taylor expanding A ca (Q ca + £) to second order in £, and using the 
third estimate in Lemma IT21 gives |M(0)| = 0oT x t\ + a~^e a e x ea if eo < 1. 

By order considerations, if we choose a — {e a e x ) s , with < s < i, and if eo << (e a e x ) 2s , then the bound 

||£|| T = 0(e a e x ) s holds for T = 0a|S|r- The bound ||£ 9 || T = 0(e a e x )^ is obtained by substituting the bound 
for into second inequality of JSJJ. We obtain the conservative estimate T\ = 0[(e a e x ) s + e t + e^]" 1 by 

substituting the bounds on ||£|| T and ||£ 9 || T into the expression for |S|t- To complete the proof, we take eo, 
e Q , and e x sufficiently small so that the smallness assumptions in all the propositions and lemmas hold. □ 

We now prove the main theorem. 

Proof of Theorem^ By our choice €q < e, there is a (maximal) time T such that the solution it in is in 
U e for time t < T . Hence decomposition 11110 with l|12|) and Proposition 1101 are valid for u and imply the 
statements of the main theorem and in particular ||£(f)||#i < Ci(e a ea;) s for times t < min{Xoj2i}. If we 
take e a ea; such that Ci(e a e.,;) s < e — 0(e a e x )i , then the above bound holds for t < T\ by maximality of the 
time Tb- d 



Appendices 

A Global Wellposedness of the bKdV 

In this appendix we prove Theorem ^ global wellposedness of the bKdV equation 

d t u = -d x {d 2 x u + u 2 +b{t,x)u) (38) 

in if 1 (R), with an appropriate condition on b. We extend the local wellposedness proof of Kenig-Ponce-Vega 
in the case of b = 0, and use an energy argument to extend local wellposedness to global wellposedness. 
Define 

\\f\\v x L< : = {j^J&tydxdt 

and similarly the L v T L q x , L^H^ and L^,W^ P norms (recall W k,p is the Sobolev space based on L P (R)). Let 
/ denote the Fourier transform of /. 

We begin with a lemma. 

Lemma 11. Let g € C(C) and g' <E L°°(C) (in distribution). If be L 2 (R) with 9 € L l (R), then 

\\\g(d x ), b]\\ L 2 {m ^ L 2 [R) <C\\g'\\ LC 
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Proof. Writting b as a Fourier integral and interchanging the integral and commutator gives 

[g(d x ), b] = — - / 6(A) \g(d x ), e ikx ] dk. (39) 



(27T) 

We require an estimate of [ff(9x)) e lfcx ] = e 4fc:E (e~ lkx g(d x )e — g(d x )). The term enclosed in brackets can 
be written as an integral: 

[g(d x ), e tkx ] = e ikx J £ (e~ ikxs g(d x )e ikxs ) ds, 
or on differentiating and computing in Fourier space, 

[g(d x ), e tkx ] = ike lkx f g'(d x + iks) ds. 
Jo 

Hence ||[5(9 X ), e ikx ] f\\ L2 < \k\ \\g'\\ Laa \\f\\ L3 for all / e L 2 (M.). Substituting this bound into proves 
the lemma. □ 

Proof of Theorem^ The proof is p power specific (/ = u p ); we only present the proof for p=2. We begin 
by proving local wellposedness. Let W be the unitary group generated by —d x ; that is, W(t) := e~ d * 1 . 
Consider the bKdV without the term bu. We rewrite this equation as a fixed point problem u — <&(«), where 
the map $ is defined by 

$ : v h-> W(t)u - / W(t - r)vd x v dr. 
Jo 

Given an initial condition u , Kenig, Ponce and Vega [2H1 proved that $ is a strict contraction on a ball 
Bx a (a) of radius a — cl{\\uq\\ Hs ) and centre v = in the space Xj, for all s > |, where 

X% := {v e C([-T,T],H S (R)) | A T ( V ) < co}, 

and 

Af» := + \\d x v\\ L ^ + IIPJ'Ml^ + (1 + T)~ p |M| w . 

Here p > |. More precisely, given ||uo||# 3 and < e < 1, there is an a and T such that uq € -Bjc s (a), 
$ : Sxj(a) -► S X| ,(ea) and Af($(«) - $(5)) < eA^(« - g). 

We now formulate Q38f) as a fixed point problem. Define 

# : u h-> - I W(t - T)d x (bv) dr. 
Jo 

Then, the bKdV is equivalent to solving u = $(w) + &(u). Thus, if $ + \& is a strict contraction on some 
ball Bx^,(a), then equation (|38|l has a solution in the same class as (|38fl with 6 = 0. To prove this we need 
estimates of Aj,(^f(v)) and Aj,(^f(v) — W(g)) for v S Bx^,(a). We will make use of the estimates 

P x \iw(t)g(x) <C||s|U, (40) 

L T L X X 

\\d x W(t)g(x)\\ L ^ L2T <C\\g\\ Ll , (41) 
\\W{t)g{x)\\ L%L ^<C{l+TY\\g\\ Hl , (42) 
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the proofs of which are given in [221 ■ 

We estimate A^(^(v)) for s — 1, beginning with || ^(f) || ^oo^i . The inequality 



W{t-T)d x {bv)dT 



< C 
+ C 



W(t-T)d x (bv)dT 



W{t-T)f) x \d x {bv)dT 



follows from the commutativity of f) x \ — —id x with W(t — r). Taking the L\ norm inside the integrals via 
Minkowski's inequality, and using that W(t — r) is unitary gives that 



\d*M\\i* l + \\PMbv)\\ Ll dr 



We increase the integration domain to [—T,T] and use Holder's inequality to obtain the estimate 



mv)\\ L¥Hh <c(\\d x (bv)\\ 



\p x \d x (bv)\\ 



The bound 



< 



\\d x W(t-r)d x (bv)\\ L „ dr 



is obtained from ||9 a; ^'(u)|| i 4 L oo by moving the derivative d x and I/°°(R) norm into the integral over r, and 
increasing the domain of integration to [—T,T]. Minkowski's inequality then implies 

\\d x *(v)\\ L *L°° < f T \\d x W(t-T)d x (bv)\\ LiL dr. 
We use that d x commutes with W and the group properties of W to rewrite this inequality: 

fT 



< 



p^W(t)W(-r)ap x \W x (bv) 



dr, 



where a is multiplication by i sgn(k) in Fourier space. The quantity W(—r)ap x \ 4 d x {bv) does not depend on 
time t. Thus, we use estimate (|40f) and that W and a preserve the L 2 (R) norm to obtain the bound 



\\d x *(v)\\ LtL „ < I p x \-*d x (bv) 

-T 



dr 



or \\d x *(v)\\ LtL „ < C \\d x (bv)\\ LtL , + \md x {bv)\ 



Again, Minkowski's inequality and the properties of W give the bound 



Wx\d x v\\ 



< 



mW(t)W(-r)d 2 x (bv) 



dr. 
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Since the same holds with integration over [— T, T] and W(— T)d%.(bv) is independent of time t, l|41l) implies 

||ftJMi~£» < \\W(-r)pMbv)\\ L}rL% < \md x (bv)\\ L , TL% . 

As above, we find that 

(l+T)-"||^)|| i2 < (l + T)-r [ T \\W(t)W(-T)d x (bv)\\ L2 L dr. 

Estimate 11421 then implies 

(l+T)-»\\*(v)\\ L2xL¥ < \\W(-r)d x (bv)\\ LlTHlx < C {\\d x {bv)\\ LlLl +\mdxibu)\\ L%L , x ) . 

Combining all estimates gives that 

A£(*(«)) < C {\\d x {bv)\\ L%L , x + \\%mbv)\\ LlLx ) , (43) 
and since ^>(v) — ^>(v) — ^(v — v) 1 

A s T (nv) - < C (\\d x (b(v - v))\\l%l\ + WPMKv - m\v T v x ) • (44) 

Thus, to prove $ + \t is a strict contraction we need estimates of the quantities appearing on the right hand 
side of (g3J) and @3J. 

Holder's inequality gives 

\\ d Abv)\\ L 2 TL 2 x < W\\ L * LJl \\v\\ L¥L 2 x + ||6|| L 2 L ™ \\d x v\\ L¥L2x , 

\\PMbv)\\ LlLl < W'\\z> t lj? Ml ¥ lI + W\\u> t l W \\9 x v\\ L¥Lx + \\P x \(bd x v)\\ L2TLx . 

We need to estimate 

\\mbd x v)\\ L 2 TLx < b\d x v\\ LlLx + \m\d x v\\ L2TLx . (45) 

The first term on the right hand side is bounded using Lemma ITT1 We obtain that 

\\m,b]d x v\\ L2TLx <c\\9\\ L2TLx \\d x v\\ L¥Ll , 

where is the L 1 norm of b' in the frequency variable. Using Holder's inequality, the second term in 

lf4"5|l is bounded as 

\\bp49xv\\ L 2 TL 2 x = \\bp4d x v\\ L 2 xL 2 T < Hftll^ioo \\Px\d x v\\ LfL 2 T . 

Combining all the estimates implies Ay(*(u)) < C||6|| X tA],(v) and A T (^(v) - *({>)) < C||6||xtA^(w - v), 
where 

||6|Ut,i := \H L 2 tW ^ + \\b\\ L 2 xL¥ + \\9\\ L 2 tL i x . 
If £ = C||6||xt,i, then the above estimates imply that VP : Bx s (a) — > Bx^(sa) and 

A^(-9(v)-V(v))<eA^(v-v). (46) 
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Thus, $ + # : B x -.(a) -> Bxj((e + £)fl) and (($ + *)(«) -($ + *) (5)) < (e + e)A^(w-w). Furthermore, if 
we take e + e < 1, then $ + ^ is a strict contraction on Bx*> (a). Invoking the fixed point theorem completes 
the proof of the local existence and uniqueness. 

Kenig, Ponce and Vega |2H] also proved that for all < T% < T, 

tik(*vo(v) - $z (v)) <C(\\v - fib|| H i +/(T 1 )(A^( U ) +A 1 Tl (v))A^ l (v-v)) , 

where f(Ti) — > as 7\ — > and <& Vo is the map associated to the fixed point problem with initial condition 
vq. The map \& is independent of initial condition; therefore, the triangle inequality and estimate (|46|l imply 

A^d^+VvM - ($ Co <C(\\v -vo\\ m + [/(2i)(A^(«)+A^(i;)) + ||6||xr]A^(i;-S)) . 

Let v be a solution to l|38[) with initial condition vq and similarly for v. Then, if T\ and ||6||xt are small 
enough, A^(w — v) < C ||i>o — ^o||#i • This proves continuity of the solution with respect to the initial 
condition and completes the proof of local wellposedness of <|38[) in iJ 1 (K). 

To extend local wellposedness to global wellposedness we require the identities 

-i poo rCfO 

d t H b (u) = - (d t b)u 2 dx and d t \\uf L2 = / b'u 2 dx (47) 

^ J— ao J -oo 

to hold for all iJ 1 (K) solutions to l|38(l . When u S H 3 (EL) both of these follow trivially by integration by 
parts. We appeal to a density argument to prove that the identities continue to hold in H 1 (M.). Let {uQ n } 
be a sequence of initial values in H 3 (M.) converging in iJ 1 (R) to Uo- Then, if (|38|) is locally wellposed in 
H 3 (M.), there are corresponding solutions u n 6 C([— T, T], H 3 (M.)) with u n (0) — ito„ and 

lim sup \\u n — u\\ H i = 0. 

n — >0 [_7^ T] x 

We have used that the time interval appearing in the local wellposedness result depends continuously only 
on the H 1 (M.) norm of the initial condition. Hence, 

1 f°° lf°° 
lim d t H b (u n ) = lim - / {d t b)u 2 n dx = - / (d t b)u 2 dx. 

n->oo rwoc I J_ OQ 2, J_ OQ 

Since dtH{u n ) — ► dtH{u) in distribution, the first identity of (|47|l holds. Similarly, 

/•OO /"OO 

lim dt \\u n \\ L 2 = lim / b'u 2 l dx= I b'u 2 dx 1 



n — >oo 



and since dt \\u n \\ L 2 — ► 9* ||u||^2 in distribution, the second identity in 1|47|) also holds. The above assumed 
local wellposedness in H 3 (R). The proof of this fact proceeds as above and one finds that (f 381) is locally 
wellposed in H 3 (R) if 

\\b\\xT,3 := \\b\\ L 2 T w^ + \\Hi? x l~ + f'h^ 

is small enough. 

We now extend the local result to a global result. The identities of (|47|) imply 



d f°° 1 f°° ( 

-j, IMIl 2 = / b'u 2 dx < e a e x \\u\\ 2 L 2 and d t H{u) = - / (d t b)u 2 dx < IM^a 
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Integrating the first by Gronwall's inequality implies ||u|| i2 < II u o|Il2 exp(e a e x t). Substituting this bound 
into the above bound on the time derivative of the Hamiltonian and integrating gives 



2 \\9xuf L 2 < \H(u Q )\ + ^- \\u a \\\ 2 exp(e a e x t) + \\uf L2 + ^ 



bu dx 



Using the bound on 1 1 1 1 ^ 2 again then gives 



1 2 / €t 1 \ > 

2 II^ w IIl2 < \H(uq)\ + I I + — + -e a ) \\uq\\ l2 exp(e a e x t) 



(48) 



This inequality implies global existence. Indeed, say there is a time T such that lim t ^T IMj^i = oo. This 
clearly contradicts ijUJ. Uniqueness follows from uniqueness of local solutions. □ 



B Proof of Lemma [3] 

Commutativity and the relation dxd^ 1 = I — ad' 1 are direct consequences of (d x + a)d~ l = I. Com- 
mutativity with S a is proved using that 

(d^ + a)- 1 :g^e- a ^ f ° g(y)e a ^ d(y - a) 



and d x + a = d x - a + a. We prove statements two and five using the above explicit formula with a = 0. 
Indeed, the inequality 

poo 

I^Vl < / \<Kx)\dx = \\4>\\ L1 



gives statement two, and since e a ^ x y>) — 1 < a\x — y\e a ( x y \ the inequality 

\^(y)\(\x\ + \y\)dydx. 

gives statement five if xcj) and ar0 are integrable. 

We prove the remaining statements in Fourier space. Let <j> be the Fourier transform of <fi. Plancherel's 
theorem implies 



\9 a 1 (p\\ L2 



(fc 2 + a 2 )- 1 |0(/c)| 2 dfc 

/CO 
(fc 2 + a 2 )- 1 (|^(fc)| 2 -|0(O)| 2 )dfc. (49) 
-oo 

The first equality immediately gives the third statement since <f> < ||^>||ri- A similar argument gives 

statement four. To prove the last statement, we concentrate on the second integral of (|49|l since the first is 
easily computed to be 

> dx' 
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When <j), xcf>, and x 2 cf> are integrable, <p and \(f>\ 2 are twice differentiable. Furthermore, since |</>| 2 is even, 
Taylor's theorem implies \$(k)\ 2 - \<j>(0)\ 2 = 0k 2 . Thus, (k 2 + a 2 )- 1 (\$(k)\ 2 - |<H0)| 2 ) is integrable for all 
a G R + and 

/oo 
(k 2 + a 2 )- 1 {\4>(k)\ 2 -\4>{Q)\ 2 )dk = 01. 
-OO 

This completes the proof. 



C Estimates of Nonlinear Remainders 



Define 

r°° 1 

N(0 :=- J F(Q ca + - F(Q ca ) - F'(Q ca )^ - -F"(Q ca )£ 2 dx 

and 

^'(0 := ~ (f(Qca + 0- /(Oca) - /'(Qca)0 • (50) 

Note that 2V'(£) = <%iV(£) under the L 2 (R) pairing. 

Lemma 12. // \\£\\ H i < 1 and / e C k (R) for some k > 3, with e L°°(M), then there are positive 
constants C\, Ci, and C3 such that 

1. \\n'{0\\l* < Ci Uf m , 

2. \\N'(0 + y"(Q ca )e\\ L2 <c 2 uf m , 

3. \N(0\ <CsUf H i. 



Proof. Taylor's remainder theorem implies 

fe-i 1 

N '® = ~ E -J (n) (Qca)C - R(Qca,Z), 

n=1 

where, since f^ G L°°(R), \R(Qca,0\ 5= C| ^ | fc . Recall that Q ca is continuous and decays exponentially to 
zero. Together with the assumption that / G C k (R), this implies f {n] {Q C a) G L°°(R) for 2 < n < fc - 1. 
Thus, after pulling out the largest constant, 

wm\v><cY l \\c\\v" 

n=2 

To obtain statementQJwe use the bound ||£ n || L 2 < C ||£||#i, which is obtained from the inequality ||£|| i0 c < 
C||^||„i and the assumption that < 1. 

Clearly, slight modification of the above proof gives items [21 and |3| For the latter we use that the 
assumptions on / imply F G C fe+1 (R) with G L°°(M). □ 
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